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Abstract. We report on the extraction of the structure functions F2 and Ax^Fj — xF!^- 
xF^ from CCFR I'^-Fe and I'^-Fe differential cross sections. The extraction is per- 
formed in a physics model independent (PMI) way. This first measurement for AxF^ , 
which is useful in testing models of heavy charm production, is higher than current 
theoretical predictions. The F2 (PMI) values measured in j/^ and /i scattering are in 
good agreement with the predictions of Next to Leading Order PDFs (using massive 
charm production schemes) , thus resolving the long-standing discrepancy between the 
two sets of data. 



Deep inelastic lepton-nucleon scattering experiments have been used to determine 
the quark distributions in the nucleon. However, the quark distributions determined 



from muon and neutrino experiments were found to be different at small values of 
X, because of a disagreement in the extracted structure functions. Here, we report 

on a measurement of differential cross sections and structure functions from CCFR 
z/^-Fe and F^-Fe data. We find that the neutrino-muon difference is resolved by 
extracting the structure functions in a physics model independent way. 

The sum of and differential cross sections for charged current interactions 
on an isoscalar target is related to the structure functions as follows: 
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Here Gp is the Fermi weak coupling constant, M is the nucleon mass, E,, is the 
incident energy, the scaling variable y = Eh/ Ey is the fractional energy transferred 
to the hadronic vertex. Eh is the final state hadronic energy, and e~2(l — |/)/(l + 
{1 — yY) is the polarization of the virtual W boson. The structure function 2xFi is 
expressed in terms of F2 by 2xFi{x, Q^) = ^2(2;, Q'^) x ^X^f^^q?) , where is the 
square of the four-momentum transfer to the nucleon, x = Q^/2MEh (the Bjorken 
scaling variable) is the fractional momentum carried by the struck quark, and 
it! = ^ is the ratio of the cross-sections of longitudinally- to transversely-polarized 
VF-bosons. The AxF^ term, which in leading order ~ 4x(s — c), is not present in 
the /i-scattering case. In addition, in a z/^ charged current interaction with s (or c) 
quarks, there is a threshold suppression originating from the production of heavy 
c quarks in the final state. For /^-scattering, there is no suppression for scattering 
from s quarks, but more suppression when scattering from c quarks since there are 
two heavy quarks (c and c) in the final state. 

In previous analyses of data, light-flavor universal physics model dependent 
(PMD) structure functions were extracted by applying a slow rescaling correction 
to correct for the charm mass suppression in the final state. In addition, the AxF^ 
term (used as input in the extraction) was calculated from a leading order charm 
production model. These resulted in a physics model dependent (PMD) structure 
functions. In the new analysis reported here, slow rescaling corrections are not 
applied, and AxF^ and F2 are extracted from two parameter fits to the data. We 
compare the values of AxFs to various charm production models. The extracted 
physics model independent (PMI) values for F2 are then compared with F2 within 
the framework of NLO models for massive charm production. 

The CCFR experiment collected data using the Fermilab Tevatron Quad- Triplet 
wide-band ly^ and beam. The raw differential cross sections per nucleon on iron 
are determined in bins of x, y, and E^, (0.01 < x < 0.65, 0.05 < y < 0.95, and 30 < 
Eij < 360. GeV). Figure 1 (a) shows typical differential cross sections at E^, = 150 
GeV. Next, the raw cross sections are corrected for electroweak radiative effects, 
the W boson propagator, and for the 5.67% non-isoscalar excess of neutrons over 
protons in iron (only important at high x). Values of AxF^ and F2 are extracted 
from the sums of the corrected t'/i-Fe and V^-Fe differential cross sections at different 
energy bins according to Eq. (1). It is challenging to fit AxF^, R, and 2xFi using 
the y distribution at a given x and because of the strong correlation between 



the AxF^ and R terms, unless the full range of y is covered by the data. Covering 
this range (especially the high y region) is hard because of the low acceptance. 

Therefore, we restrict the analysis to two parameter fits. Our strategy is to fit 
AxFs and 2xFi (or equivalently F2) for x < 0.1 where the AxF^ contribution is 
relatively large, while constraining R using the R^^orid QCD inspired empirical fit 
to all available R from electron- and /^.-scattering data. The R'^orid ^^^^ 
good agreement with NMC i?'^ data at low x, and with the most recent NNLO 
QCD calculations (including target mass effects) of R by Bodek and Yang 

For X < 0.1, J? in neutrino scattering is expected to be somewhat larger than 
R for muon scattering because of the production of massive charm quarks in the 
final state. A correction for this difference is applied to R'^^orid using a leading 
order slow rescaling model to obtain an effective R for neutrino scattering, -Re//- 

The difference between R'^^orid R'^^f is used as a systematic error. Because 
of the positive correlation between R and AxF^, the extracted values of F2 are 
rather insensitive to the input R. If a large input R is used, a larger value of 
xF^ is extracted from the y distribution, thus yielding the same value of F2. In 
contrast, the extracted values of AxFs are sensitive to the assumed value of R, 
which is reflected in a larger systematic error. The values of AxF^ are sensitive to 
the energy dependence of the neutrino flux (~ y dependence), but are insensitive 
to the absolute normalization. The uncertainty on the fiux shape is estimated by 
using the constraint that F2 and xF^ should be fiat over y (or E^) for each x and 
g2 bin. 

Because of the limited statistics, we use large bins in in the extraction of 
AxF^ with bin centering corrections from the NLO Thorne & Roberts Variable 
Flavor Scheme (TR-VFS) calculation with the MRST PDFs. Figure 1 (b) shows 
the extracted values of AxF^ as a function of x, including both statistical and 
systematic errors, compared to various theoretical methods for modeling heavy 
charm productions within a QCD framework. The three-flavor Fixed Flavor Scheme 
(FFS) assumes that there is no intrinsic charm in the nucleon, and all scattering 
from c quarks occurs via the gluon-fusion diagram. The concept behind the Variable 
Flavor Scheme (VFS) proposed by ACOT is that at low scale, /x, one uses the three- 
flavor FFS scheme, and above some scale, one changes to a four-flavor calculation 
and an intrinsic charm sea (which is evolved from zero) is introduced. The concept 
in the RT-VFS scheme is that it starts with the three-flavor FFS scheme at a low 
scale, becomes the four-flavor VFS scheme at high scale, and interpolates smoothly 
between the two regions. Shown are the predictions from the TR-VFS scheme (as 
corrected after DIS-2000 and implemented with MRST PDFs), with their suggested 
scale fi = Q, and the predictions of the other two NLO calculations, ACOT- VFS 
(implemented with CTEQ4HQ and the recent ACOT suggested scale /i = rric for 
Q < rric, and //^ = rric^ + cQ'^{l — nic^/Q'^Y for Q < rric with c = 0.5 and n — 2), 
and the FFS (implemented with the GRV94 PDFs and GRV94 recommended scale 
/i = 2mc). Also shown are the predictions from AxF^ ~ 4Ks{x, Q"^) from a leading 
order model (LO(CCFR)) Buras-Gaemers type fit to the CCFR dimuon data (here 



K is a slow rescaling correction). Figure 1 (b) (right) also shows the sensitivity to 
the choice of scale. The data do not favor the AC0T-VFS(CTEQ4HQ) predictions 
if implemented with an earlier suggested scale of fi = 2Ptmax- With reasonable 
choices of scale, all the theoretical models yield similar results. However, at low 
our AxF^ data are higher than all the theortical models. The difference between 
data and theory may be due to an underestimate of the strange sea (or gluon 
distribution) at low Q^, or from missing NNLO terms. 

As discussed above, values of F2 (PMl) for a; < 0.1 are extracted from two 
parameter fits to the y distributions. In the x > 0.1 region, the contribution from 
AxFs is small and the extracted values of F2 are insensitive to AxF-^. Therefore, 
we extract values of F2 with an input value of R and with AxF^ constrained to 
the TR-VFS(MRST) predictions. As in the case of the two parameter fits for 
X < 0.1, no corrections for slow rescaling are applied. Fig. 2 (a) shows our F2 
(PMI) measurements divided by the predictions from the TR-VFS(MRST) theory. 
Also shown are F2 and F2 from the NMC divided by the theory predictions. In 
the calculation of the QCD TR-VFS(MRST) predictions, we have also included 
corrections for nuclear effects, target mass and higher twist corrections at low values 
of Q^. As seen in Fig. 2, both the CCFR and NMC structure functions are in good 
agreement with the TR-VFS(MRST) predictions, and therefore in good agreement 
with each other. A comparison using the AC0T-VFS(CTEQ4HQ) predictions 
yields similar results. 

In the previous analysis of the CCFR data, the extracted values of F2 (PMD) 
at the lowest x — 0.015 and Q'^ bin were up to 20% higher than both the NMC 
data and the predictions of the light-flavor MRSR2 PDFs. (see flgure 2 (b) ). 
About half of the difference originates from having used a leading order model for 
AxF^ versus using our new measurement. The other half originates from having 
used the leading order slow rescaling corrections, instead of using a NLO massive 
charm production model, and from improved modeling of the low PDFs (which 
changes the radiative corrections and the overall absolute normalization to the total 
neutrino cross sections). 

In conclusion, the F2 (PMI) values measured in neutrino-iron and muon- 
deuterium scattering show good agreement with with the predictions of Next to 
Leading Order PDFs (using massive charm production schemes), thus resolving 
the long-standing discrepancy between the two sets of data. The flrst measure- 
ments of AxF^ are higher than current theoretical predictions. 
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FIGURE 1, (a) Typical raw differential cross sections at Eiy = 150 GeV (both statistical and 
systematic errors are included), (b) Axi^3 data as a function of x compared with various schemes 
for massive charm production: RT-VFS(MRST), AC0T-VFS(CTEQ4HQ), FFS(GRV94), and 
LO(CCFR), a leading order model with a slow rescaling correction (left): Also shown is the 
sensitivity of the theoretical calculations to the choice of scale (right). 
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FIGURE 2. (a) Left side: The ratio (data/theory) of the FJ^' (PMI) data divided by the predic- 
tions of TR-VFS (MRST) (with target mass and higher twist corrections). Both statistical and 
systematic errors are included. Also shown are the ratios of the F2 (NMC) and Fl (SLAC) to the 
TR-VFS(MRST) predictions, (b) Right side: The ratio (data/theory) of the previous (PMD) 
data (and also F^ (NMC) and F| (SLAC)) divided by the predictions of the MRSR2 light-flavor 
PDFs (with target mass and higher twist corrections). 



